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We calculate the equation of state of a gas of strings at high density in a large toroidal universe, 
and use it to determine the cosmological evolution of background metric and dilaton fields in the 
entire large radius Hagedorn regime, (InS)^/'* < 7i < 5^/'* (with 5* the total entropy). The 
pressure in this regime is not vanishing but of 0(1), while the equation of state is proportional 
to volume, which makes our solutions significantly different from previously published approximate 
solutions. For example, we are able to calculate the duration of the high-density "Hagedorn" phase, 
which increases exponentially with increasing entropy, S. We go on to discuss the difficulties of the 
scenario, quantifying the problems of establishing thermal equilibrium and producing a large but 
not too weakly-coupled universe. 



I. INTRODUCTION 

Strings may have played an important role in the 
very early universe due to the potential dominance of 
stringy effects at high energies. It is important therefore 
to understand the dynamics of an ensemble of strings in 
what seems to be a rather natural background, namely 
dilaton gravity, the stringy version of Einstein gravity. 
This approach to string cosmology started with the work 
of Brandenberger and Vafa (BV) [l| who laid down the 
main conceptual framework and subsequently Tseytlin 
and Vafa (TV) [3] who introduced dilaton gravity into 
the picture, hence making a first step towards realising 
a string cosmology at finite temperature. This frame- 
work incorporates all the degrees of freedom of pertur- 
bative string theory, namely oscillator, momentum and 
winding modes, which at high densities produce a state 
with unusual thermodynamic properties such as a lim- 
iting temperature, known as the Hagedorn phase. An 
early Hagedorn phase may turn out to have crucial con- 
sequences for the early universe, the initial singularity 
and the dimensionality of spacetime (see [l| for further 
details). 

Even though this scenario has been studied and ex- 
tended in a number of different ways (see eg [1, i, S 
B 0, B i, 0, [H E E H), a complete string theory 
analysis has yet to be made. Furthermore, a viable string 
cosmology should also incorporate the attractive features 
of the standard cosmological model, such as the absence 
of a flatness and horizon problem and solutions to the 
entropy and size problems. 

We shall here discuss the scenario introduced by TV 
which was inspired by the conceptual framework of 
BV We shall concentrate on large radius evolution, 

< i? < S'^/'^ (with S the total entropy), where the 
dominant contribution to the pressure is from momen- 
tum modes (roughly speaking the centre of mass motion 
of small loops) and hence not negligibly small. In par- 



ticular, we shall be tracing the evolution of the entire 
large radius Hagedorn phase up until R ~ S^/'^ beyond 
which the universe is radiation dominated. The large ra- 
dius regime had been poorly studied to date and only 
few analytic results had been obtained. Furthermore, 
analysing the large radius dynamics is crucial in order to 
understand various transitions that occur and to make 
contact with late time evolution (by which we mean ra- 
diation domination onwards). Finally, a large radius uni- 
verse cannot be avoided in the dilaton gravity realisation 
of String Gas cosmology, making the reasonable assump- 
tion that the massless string degrees of freedom which 
take over at the end of the Hagedorn phase are identified 
with the particle content of the Standard Cosmology. A 
simple calculation then shows that the compactification 
radius must be at least 10~^ mm in order to expand to 
at least the Hubble radius today. 

It has been proposed that the string gas initial con- 
ditions can produce scale- free cosmological perturbations 
15 , [l^ , although the proposal was strongly criticised in 
17| . In order to rescue the scenario more recent vari- 
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ants depart from the original simplicity of the scenario 
by fixing the dilaton and introducing extra terms and 
fields [l3, Ell • Another variant supposes that before 
the perturbative Hagedorn phase of the original scenario 
there was a strongly-coupled Hagedorn phase with con- 
stant dilaton [23|. Our analysis applies equally to the 
perturbative Hagedorn phase of this scenario - this is 
independent of whether the dilaton is initially strongly 
coupled provided the low energy effective action is valid 
when the universe is much larger than string scale. 

In section |TT] we give a brief overview of the deriva- 
tion of the equations of motion of string cosmology to 
set our notation conventions. In section [TTTl we reformu- 
late the equations of motion so as to make them analyt- 
ically solvable and discuss their general structure while 
making explicit the constraints that appear. We work in 
the string frame but the Einstein frame evolution is also 
briefiy discussed, especially when making contact with 
late time evolution and observational constraints. In sec- 
tion IIVI we derive the equation of state parameter which 
(to leading order) turns out to be proportional to the 
inverse entropy density. We then go on to derive the 
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number density associated with winding and momentum 
modes of string, which correspond roughly speaking to 
long and short string loops respectively. We notice that 
at radii larger than ~ (In 5')"'^/'* the contribution of long 
string to the number density is negligible compared to 
that of small loops while at radii greater than ~ S'^/'' 
there is not enough energy to excite massive modes of 
string and only zero modes are present, i.e. radiation. 
In section |V] we solve the equations of motion analyti- 
cally with evolving equation of state for the entire large 
radius Hagedorn regime and compare our results with 
the corresponding numerical solutions. We find that ne- 
glecting the evolution of the equation of state cannot be 
justified in the large radius Hagedorn regime and em- 
phasise that the pressure is not negligible but is rather 
0(1) in string units and approximately constant. The 
problems of this scenario are discussed in section IVIl In 
particular, the issues of Jeans instabilities and maintain- 
ing thermal equilibrium are examined and we improve on 
previous estimates H, [U by showing that the universe 
could not have been in thermal equilibrium at any stage 
of the evolution in the large radius Hagedorn region if 
string is weakly coupled, i? <C 1. We also discuss the size 
problem, and its relation to the weak-coupling problem 
in String Gas cosmology. 



II. GENERAL FORMALISM 

The critical tree level dilaton-gravity low energy ef- 
fective action plus stringy matter contributions reads [s^l 
(see eg [13) 

r^xR 

X [Rw + 4G^^Va<^Vb0 +...]+ 5m, (1) 

where (f> is the dilaton, Riq is the 10-dimensional string 
frame Ricci scalar and Gab with A, i? = 0, • • • , 9 is the 
string frame metric tensor [ssj . The first two terms rep- 
resent the lowest order (in both a' and gs) contribution 
to the bosonic NS-NS sector and are present in all 5 su- 
perstring low energy effective actions [s^l- The string 
gas contribution is contained in the matter action, Sm, 
the form of which is derived from the string spectrum on 
a torus at finite temperature. We are thus coupling a 
gas of strings to a dilaton-gravity background. In writ- 
ing down the above action we have assumed that the 
string coupling is small, ~ e^"^ ^ 1, that the curva- 
ture scalar Rioa' <C 1 and that field moduli evolve slowly 
(adiabaticity). There are branches of solutions for which 
these conditions are attractors as we shall demonstrate 
when the matter contribution, that is quantified by the 
equation of state parameter, w, varies slowly with size 
(see also 0). The quantity that we cannot neglect is the 
dilaton because this is required for the T-duality invari- 
ancc of the theory. The equations of motion then follow 



from varying S = S[Gab, (p]' 

R^B + 2V'^Vb<^ = ^ae'^T^, (2) 

i? - 4(V0)2 4V20 = 0. (3) 

We shall henceforth work in units where k\q = 1/2. Let 
us concentrate on homogeneous evolution of metric and 
on homogeneous and isotropic evolution of dilaton 

ds2 = -dt2+^e2A,(t)^^2^ = 0(0. (4) 
1=1 

The energy-momentum tensor is effectively, as we shall 
demonstrate in section HVl that of an ideal fluid but with 
variable equation of state, T"^^ = diag(— p, pi, • • • jPg). 
The equations of motion ^ [3]) then reduce to [1] (see 
also |23|): 

-^Af + vi^=e^+^.^V, (5) 
1=1 

A.-^A, = ie'^+^.^^K, (6) 

i=l 

where (p is the re-scaled dilaton 

9 

(p=20-^A,. (8) 

i=l 

We are neglecting the effect of interactions and so we are 
studying the evolution of the universe after the hierarchy 
of scales has been created due to the BV mechanism. We 
are interested therefore in the case when d dimensions are 
large, equal and free to expand and are assuming 9 — c? 
dimensions are static and at string scale (the case d = 3 
being of interest for the BV scenario). To see that this 
factorisation is consistent with the equations of motion 
note that at the string scale the pressure is negligibly 
small [ll,[l,01. Then, writing 

r A(i) for * = {!,..., 4, 
^'^^>~\^x{t) fori = {d+l,...,9}, 

it follows from ^ that /i(t) = const is a particular so- 
lution when Pi\i>d+i = 0. However, the pressure, pi, is 
vanishing only at string scale and so it is consistent to 
take 

In making this deduction we have assumed that the equa- 
tions of motion presented above hold at string scale and 
so this may not be a good assumption because non- 
perturbative effects are expected to become important 
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at the string scale. The density and pressure of the 
large dimensions relate to the enery via p{X) = E/V 
and p{X) ~ —dE/dV respectively with V = e'''^. Only 
two of the equations (071) are independent (([5]) being a 
constraint). These equations have been studied 0, 0| 
in some detail for p ~ (string scale Hagedorn region) 
and p = p/d (massless momentum mode or radiation 
domination) and so we shall concentrate on large radius 
Hagedorn regime where the pressure, p, is approximately 
constant and (as we shall show in section|V]) the variation 
of the equation of state cannot be neglected. 

In the next section we recast the equations of motion 
into a form that makes them easier to solve analytically 
and determine their underlying structure. 



III. EQUATIONS OF MOTION IN STRING 
AND EINSTEIN FRAME 

A. String Frame 

Let us identify w{\) with the equation of state pa- 
rameter p = w{X)p. In section lTVl we shall derive the form 
of w(A) from the microcanonical ensemble. This defini- 
tion enables us to reduce the equations of motion to first 
order uncoupled form by introducing the two following 
variables 



Y{\) = dip/dX and H{X) = X, 



(10) 



where H is evidently the (string frame) Hubble parame- 
ter. With these definitions the equations ([5][7]) on account 
of ([9]) lead to the following two independent equations of 



motion 



dY 

Ix 



(1 + wY) (y2 _ d) , 



(11) 



(12) 



The constraint ([5]) will serve to set the initial conditions 
and takes the form 



.20 



Y^-d' 



(13) 



It thus becomes clear that the dilaton. its rate of change 
and the number of dimensions determine the evolving 
gravitational coupling. The structure of pip is pictorially 
depicted in Fig [1] while that of is shown in Fig [51 

The evolution of the dilaton follows immediately 
from (El) and (HOl): 



(14) 



To determine the constraints on Y note that from the 
constraint equation we learn that positivity of the energy 
implies 



p > ^ y2 > d. 



(15) 



The equations of motion tell us that this condition is at- 
tractive. Furthermore, in order for the action ([1]) to be 
justified we require the branch of solutions that docs not 
allow runaway increasing dilaton in an expanding uni- 
verse: Y < —\/d. In fact the decreasing dilaton branch 
of solutions satisfies, 



< ^ 



Y <-d \iH > 0, 

Y > -d if i^ < 0, 



(16) 



which follows from the definitions (fTUll of H and Y . The 
condition for an increasing dilaton, (/> > 0, can be reached 
by reversing the inequalities involving H. Notice that an 
increasing dilaton is possible in an expanding universe 
when Y > —d but when Y < —^/d it reaches a maximum 
and then decreases as the Y = —d line is crossed. 

The equations are invariant under T-duality: A > 
—A, as T : y 1-^ —Y and T : w ^ —w. Hence, we shall 
consider only the w > branch as the dual branch w < 
can be reached with a T-duality transformation. Wc also 
conclude from (fTTj) that 



Fad (A) 



wiX)- 



(17) 



is an attractive solution when w(A) varies sufficiently 
slowly with increasing A. We shall refer to (fTTl as the 
adiabatic solution. By sufficiently slowly we mean that 
w~^dw/dX <C 1 should be satisfied. In section |Vl we 
show that this condition is not satisfied in the large ra- 
dius Hagedorn regime (in fact w~^dw/dX ^ d/w ^ 1) 
and so the adiabatic approximation is not a good approx- 
imation. It will prove useful however in that it enables 
us to obtain analytical solutions by expanding Y about 
— 1/w. See also (fTBl) and the dash-dot line in the fiow 
diagram Fig [H from which the attractor solution can be 
read off. Notice also that this solution corresponds to the 
minimum of the dash-dot line in Fig [51 where it is seen 
that close to the adiabatic solution the Hubble parameter 
decreases in an expanding universe. 

At radiation domination (jj approaches a constant (as 
for radiation w = 1/d, see also (fT4l) ) [ssj . 0' — > 0, and 
hence Y ^ —d (the dash-dash line in Fig |ll2p . 



B. Einstein Frame 

The string frame metric is conformally related to 
related to the Einstein metric via 



G 



flu 



exp 



1 



G 



(18) 



as can be verified by direct substitution of the new metric 
into the action ([T]) . Wc shall throughout denote Einstein 
frame quantities with a tilde. It follows that the Einstein 
and string frame Hubble parameters arc related by 



H 



24> 

gd-1 



-Y - 1 



(19) 
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FIG. 1: A flow diagram for dY/dX. The dash-dot curve 
represents the dependence of dY/dX on Y for some value 
< w{X) < 1/d, while the dashed line represents this relation 
for w = 1 /d. This corresponds to the asymptotic form of Y' vs 
Y in an expanding universe where Y < —d, see (|15I16|) . The 
solid curve corresponds to the maximum value of w allowed 
by the dominant energy condition, < 1. At Y = ±\/d 
(represented by the squares) the value of Y' is independent 
of w. We see that the solution Y = —1/w (represented by a 
circle) is attractive for all initial Y < —y/d (provided w varies 
slowly with size) and tends to — oo close to the self dual point 
where w =^ 0. 



from which, on account of piH2p . the rate of change of 
H with respect to A follows 



= -d+-(d - 1)(1 - w{\)) ^^ 1 , 



(20) 



where A = A — 2(f) /{d — 1). Notice that we have writ- 
ten (|20p in terms of string frame quantities Y {X),w{\). 
We shall mainly concentrate on string frame calculations 
and we have written down p8ll9|20p to gain an idea of 
how the corresponding Einstein frame quantities evolve. 
From (fin)) and pB]) it follows that in the decreasing dila- 
ton branch of solutions expansion in the string frame ne- 
cessitates expansion in the Einstein frame. 

In the next section we derive the thermodynamic 
observables of interest starting from the microcanoni- 
cal density of states. We concentrate on large compact 
spaces and derive the dependence of z«(a;) on x (~ R'^/S, 
the inverse entropy density) p5|) that will enable us to 
solve the equations of motion in the regions of interest. 



FIG. 2: A sketch of (InH)' vs Y (equation The dash- 

dot curve represents this relation for some value < w{X) < 
1/d, while the dashed line represents this relation for w — 1/d 
(the asymptotic form of (ln_H")' vs Y in an expanding universe 
where Y < —d (see (|15I16|) )). The solid curve corresponds to 
the maximum value of w allowed by the dominant energy con- 
dition, |«;| < 1. At y = ±\/d (represented by the squares) the 
value of (In H)' is independent of w. From Fig[T]we know that 
the solution Y = —1/w (represented by a circle at the mini- 
mum of the solid curve) is attractive for all initial Y < —\fd 
provided w varies slowly with size. This minimum propa- 
gates from Y — > — oo, when ixj ~ 0, through to K — > — d, 
when w = 1/d, in which case the dashed and solid curves 
coincide. Whether w > 1/d is possible is not clear at present 
but the leading terms of the density of states seem to imply 
that w > 1/d is possible. Note that when the density is below 
the Hagedorn density radiation is expected to dominate the 
spectrum and so to = 1/d eventually. 



gion for both small and large spaces, compact and non- 
compact topologies. Here we merely present the results 
for the case of a large toroidal background in the vicinity 
of the Hagedorn temperature, where the string density is 
higher than 1 in string units. 

The microcanonical approach starts from the defi- 
nition of the microcanonical density of states ^l{E) = 
J2a ~ Ea), with the sum taken over all states a of 
total energy Ea ■ We shall concentrate on the case where 
all 9 spatial dimensions are toroidally compactified with 
9 — d dimensions at string scale and the remaining d di- 
mensions of equal size i? 3> 1 . 

The leading two terms of the total density of states 
for a large and compact stringy universe when E ^ piR'^ 
turns out [1^ to be given by. 



IV. SUPERSTRING THERMODYNAMICS 

Thermodynamics of non-interacting superstring 
gases was studied extensively in the late I980's but the 
work that is most relevant for the present paper is found 
in raj|Jp,[2i,[23,[Ii,[li,[33,[3l|,l33 and more recently 
in [33l. |34|. Focus has mainly been on the Hagedorn re- 



[{(3o-(3,)Er-' 

{2d-iy. 



-(Po-Pi){E-piR'') _ ^ 



(21) 



where i? ^ I and d > 2. The constants Cq, bo and pi are 
0(1) and pi is identified with Hagedorn energy density. 
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The quantities /3„ (the cases n = 0, 1 being of immediate 
interest) are defined to be the inverse temperatures at 
which the analytically continued (single or multi-string) 
partition function is singular. In particular, the leading 
singularities for d large equal radii were found to be lo- 
cated at [23, Hi 



large volume string gas we find that the pressure 



/3„=(2^V)'/'x ^ ^ 

... (22) 

for positive or vanishing integer n subject to n/2i?2 < 1 
and R = Rj^fa' (note that /3o > /^n^^o)- The inverse 
temperature /3 = /3o is by definition the inverse Hage- 
dorn temperature. Therefore, Po — Pi — /3o/(4\/2^^) 
when 1/2R^ <C 1. Note that the second term in the 
brackets of pT|) differs from [2^, and a more detailed 
derivation is given in [11] . (The right hand side of equa- 
tion (3.15) in [2^ should have an additional factor of 

The leading behaviour fl ~ e^°^ is well known in 
string theory and is due to the exponentially rising de- 
generacy of states associated with high energy oscillator 
modes of string (see e.g. [l^). An important differ- 
ence between (|2ip and the corresponding expression for 
all radii at string scale is the factor gqR'^ that appears 
in the exponent. This is absent when all dimensions are 
of stririg scale and is due to the momentum modes of 
string [32| . the corresponding winding modes are sub- 
leading and their contribution can be understood by a 
T-duality transformation on f2. In particular, the corre- 
sponding expression for the density of states for an ideal 
gas of light point particles [l^l is exp(p'*/(''+i)i?''). 
The proportionality factor in the case of particles is the 
density, p, and this varies with size. The reason as to 
why the coefficient of i?'^ is constant in the string calcu- 
lation (j21[) is that the momentum modes density is con- 
stant above Hagedorn energy densities (see also U for a 
classical simulation that supports this view). It is hence 

clear that the term e°°^ appears due to the presence of 
momentum modes, roughly speaking the centre of mass 
motion of short loops of string. 



n \ 



1/2 



2R^J 



1/2 



A. Equation of state 

From the entropy S* = InO (up to an overall con- 
stant), all thermodynamic observables follow. From the 
assumption of adiabaticity, S = hiQ{E, R) = const, with 
il{E,R) given by (HU, we find that total energy of the 
string gas is 



E{x) = Po^S{l - .T - 0{x/R)), 



(23) 



where we have found it convenient to define x = a^R'^/S, 
which is essentially the inverse entropy density. This ex- 
pression is valid for i? ^ 1 and .t ^ 1. 

The pressure of the string gas system is given by 
p = p-^{dS/dV)E, with V = R'^. In particular, for the 



p{x) = QoPq ^ + 



d-l bp _ 2p [3o - pi 
d R d Po 



(50 



(24) 



where 



(50 



[{Po - Pi)E] 



^-iPo-MiE-piR") 



(2d - 1)! 



This result strictly applies in the region p ^ pi. Note 
however that the second term in the expression for the 
pressure is a correction that arises due to a higher order 
term in the leading exponential of (PT|) while the third 
term is associated to (50. In particular, there are in prin- 
ciple other corrections that are subleading relative to the 
second term but larger than the third term which also 
come from the leading exponential in (|2ip . The reason we 
have written the expression for p{x) in the above form is 
because the corrections from the higher order terms in the 
exponential decrease with increasing radius whereas (50 
increases. Therefore, a careful treatment of the sublead- 
ing contributions is needed, especially close to the Hage- 
dorn density, and we hope to pursue this in our follow-up 
paper [s^ . For the current paper the important term is 
the leading contribution while the sub-leading corrections 
have only been written down to obtain an order of magni- 
tude estimate of the region of validity of the leading term. 
It is important to emphasise that the pressure does not 
vanish in this large radius Hagedorn regime (even though 
p/p ^ x) but rather is of order qoPq^ . We find that the 
equation of state parameter, w = p/p, is correspondingly 
given by 



w{x) 



1 



1 + (bR-^) 



1 



2 po~Pi 
d Po 



(50 



(25) 

where b is an 0(1) constant. This expression holds for a 
non-interacting string gas and breaks down at a critical 
volume Vc ^ S (or Xc ~ 1) in which case (50 ~ 1. In 
what follows we shall only consider the leading term. 

For very large volumes the universe will be radiation 
dominated and w ~ 1/d. The equation of state in the 
Hagedorn/radiation transition region a; ~ 1 has not been 
determined to date, the calculation of which boils down 
to determining a sum of contour integrals of the analyt- 
ically continued multi-string partition function around 
each of the singularities /3„. Some progress is made in 
this direction in [35| . Furthermore, we have mentioned 
that the pressure close to the string scale, i? ~ 1 (in 
string units), vanishes for the string gas whereas we have 
shown in ((24)) that for large radii it is of order 1. The 
dependence of the string gas pressure on size between 
these two extreme cases is also puzzling and has not been 
determined to date. However, we would rather not say 
much about the small radius region where the validity 
of the free string approach is questionable. In particu- 
lar, one would expect strong coupling effects to become 
important. 
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B. String number density 

One can also estimate the number density of small 
loops and long winding strings. The total number of 
strings in an ensemble with individual energies in the 
range e ^ e + de is classically given by [l^l 



V{e;E) de 



1 



f{e)n{E-e) de. 



(26) 



/(e) is the single string density of states, /(e) = J^a ^{^^ 
en) and e^ is the energy of the single string state a. In 
[26| it is shown that (for single string energies greater 



than some lower cut-off eo 



-1/2 



(27) 



where the multiplicity of states, g„, is defined to be the 
total number of possible configurations for which n ~ 

SiLi "T-l foi' ™i £ ^- For example, go = 1, 5i = 2d, 
52 = 2^d{d - l)/2!, 53 = 2^d{d - l){d ~ 2)/3! and so 
on, whereas for very large n we can make a continuum 
approximation and find 



2Tr^ 
r{d/2) 



,d/2 



(28) 



This corresponds to the number of states contained in a 
spherical shell of thickness 1 and radius n^/^ [s^. For the 
rest of this section we find it convenient to set a' = 1. 

Equation (|27p can then be evaluated for large and 
small single string energies leading to [1^ 



fie) 



cd/2+l ' 



eo < e < i?^ 



e>i?^ 



(29) 
(30) 



where eo is some cut-off (~ 1) below which our expression 
for density of states breaks down as mentioned above. 
The first of these is calculated using the continuum ap- 
proximation for the density of states, gn 



where n 

is large. This approximation can be made because, as 
can be seen from (j^ . the radius-dependent singularities 
Pn>o all ap proach Pn (from below) as the radius R in- 
creases |25| thus making the monotonic map /3„ : Z — > K 
approximately continuous and slowly varying. Therefore, 
all factors of the form e^">°'^/e will be comparable for 
arbitrary n provided R^ > 2n, and so given that gn in- 
creases with increasing n we expect the large n region 
of gn to dominate the sum in (j27p where the continuum 
approximation (j28p can be made. The sum in (|27p is 
therefore to be replaced by an integral which can be per- 
formed explicitly to yield ((^ . The second of these ((501) 
is calculated in the small radius limit. Here the lead- 
ing (Hagedorn) singularity is far from any sub-leading 
singularities /3„>o so that terms such as e^">°'^/e are ex- 
ponentially suppressed compared to e^"'^/e. Hence, only 



the leading {n = 0) term of the sum in ((27)) remains for 
large single string energies ((50]). 

The high energy strings with e ^ R^ correspond 
to long winding strings which can accommodate a larger 
number of oscillators compared to small loops which in 
turn correspond to low energies e <^ R^. Note further- 
more that one can interpret states with e ~ i?^ as random 
walks of step length 1, whose size is of order R. 

We can estimate the total number of strings, N{E), 
in the ensemble from 



N{E) 



1 



^(E) JO 

defie)e-^'>' 



def{e)n{E~e) 



(31) 



Let us then write this integral as a sum of two contribu- 
tions on account of ((^5)) and ((50)) to find that 



R'^ + \n{E/R^ 



(32) 



Therefore, from ((32p and ((23)) we find that the number 
density of strings for a space with d large compact di- 
mensions will be 



n{R) ^l + R-'^\nS. 



(33) 



Roughly speaking, the first term counts the number den- 
sity of small loops (coming from the e <C i?^ term) while 
the second term counts the long string contribution (com- 
ing from the e '3> R^ term) (see |2a]). There is hence 
about one loop per unit volume and we see that long 
strings are energetically favourable at radii smaller than 
~ (InS)^^'^. For radii larger than ^ (InS)^^'^ small loops 
will dominate the number density and the contribution 
from long string is expected to be small. Furthermore, 
at a critical radius Rr S^^'^ there is not enough en- 
ergy to excite the massive states however and massless 
momentum modes (radiation) will come to dominate [l| 
the ensemble. This radius, Rr, is reached when the total 
energy density is of order one, oy E ^ PiR"^, and it is 
seen from ()21() that this is also where our calculations 
break down; in fact higher order terms in the expansion 
of Q also become important at this density as one would 
expect [35j . Notice that the number density of windings, 
n.w, will be completely negligible at this critical radius 
Rr ^ S^''^: 



(ln^)/5< 1. 



(34) 



We have thus found that there are three distinct 
regions in the large radius regime. For radii R in the 
range 



K i?< (InS-) 



l/d 
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there is a large number of long winding string compared 
to small loops. For radii in the range 

(In 5)1/'^ < i?< 5!/^ 

the large radius Hagedorn region that is most relevant for 
the rest of this paper, there is expected to be only a small 
number of long string compared to small loops which 
dominate the number density. Note that this transition 
from long to short string does not affect the equations of 
state. Finally, for radii R in the range 

S^/"^ < R, 

i.e. outside the large radius Hagedorn region, there is not 
enough energy to excite massive modes and the spectrum 
is dominated by the ground state which is comprised of 
massless modes, i.e. radiation. This latter transition 
does affect the equation of state. 



where " . . . " denote terms of higher order in w. Note that 
e{w) is not necessarily small and that \w\ < 1. From 
we see that the differential equation for the string frame 
Hubble parameter H can correspondingly be written in 
terms of e(w) as follows 



dw 2(1 + 



1 



(38) 



The first thing to notice is that there are two opposing 
terms in p7p . If e/w 2d then there is a runaway 
solution so that the adiabatic solution e = is repulsive. 
On the other hand, when e/w ^ 2d the e = solution 
becomes attractive. Independently of which of these two 
initial conditions the universe chooses the attractor of 
P7p has approximately the following form 



EVOLUTION IN STRING FRAME 



Vl + 8dw - 1 



(39) 



In this section we shall present solutions for the Hub- 
ble parameter and the dilaton as a function of size for the 
equation of state parameter ((25)) that was derived in the 
previous section. Keeping only the leading contribution 
we have 



w(x) 



(35) 



Recall that x = a^R'^/S <JC 1 is the inverse entropy den- 
sity multiplied by an 0(1) constant ag > 0. 

As a first approximation let us consider the adiabatic 
attractor solution ()17p (in an expanding universe with 
decreasing dilaton) 



It becomes clear that we cannot identify w^d (x) with (pS)) 
because the condition that we trust the adiabatic approx- 
imation, w^'^\dw/dX\ <C 1, is violated for the equation of 
state parameter ([35)) . for which w^'^dw/dX ~ d/w 3> 1. 
Furthermore, this is true for the entire range oiw allowed 
by the dominant energy condition. Let us then expand 
about the adiabatic approximation in the following man- 
ner: 



Y(x) 



1 



w{x) 



{l + e{x)). 



(36) 



On account of the equation of motion of Y{\) (fTT)) . we 
obtain the equation of motion of e{x). We shall find that 
e increases with increasing w and shall hence refer to the 
solution (|36p as the non-adiabatic approximation. We 
find it convenient to take as our independent variable 
the equation of state parameter w; on account of (|55)) . 
X = w/{l + w). Then, the equation of motion for e{w) is 



de{w) 
dw 



1 



1 - 



1 e(l 



2dw{l + w) 



(37) 



Furthermore, it follows from (|55)) that when e ^ 1 + 
0(1/7^) the Hubble parameter increases with increasing 
size. This, on account of (P5)) . will be the case when 
w > 1 /d. In other words, if w overshoots the momentum 



mode value, the string frame Hubble parameter will 
increase. 

We have shown that the true evolution departs from 
the adiabatic approximation, Y ~ —l/w, by a factor of 
approximately (Vl + 8dw + l)/2. We can however derive 
the evolution of the dilaton from the definitions (|55)) and 
(HD) for the solution (|55)) and we find that 



Cj){w) ~ 0(«7o) + 



VI + 8dw + 1 
4dw 



In 



VI + 8dw + l 
Vl + 8dw 



1 



9.' AH 



1 



2d 



arctan 



l + 8dw 



1 



fiiwo), (40) 



where /i(u'o) is an integration constant such that 
4'{'>^)\wo — '^(w'o)- This equation tells us that the dila- 
ton decreases rapidly and monotonically with increasing 
equation of state parameter w. In the limit w I the 
dilaton decreases according to 



(w) - 0(wo) 



1 



1 



1 



2d \wo 



(41) 



where " . . . " denote higher order terms in w. The solution 
(|40)) is plotted alongside the numerical solution (by which 
we mean the numerical solution to equation ()lip when 
sourced by (|35p ) in Fig [3] We can likewise determine the 
evolution of the string frame Hubble parameter and from 
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FIG. 3; A plot of the dilaton, (f>{x), versus x for wo — 0.1. 
The dash-dash Une corresponds to the full numerical solution 
while the solid line corresponds to the analytic solution (|40[). 




FIG. 4: A plot of the Hubble parameter, H{x), versus x for 
Wo = 0.1. The dash-dash line corresponds to the numerical 
solution while the solid line corresponds to the analytic solu- 
tion (|42|) . Notice that there is a minimum H{x) when w ~ 1/d 
so that H actually slightly increases beyond this value until 
X ~ 1/2 where to ~ 1. 



dSHl), (I2S1) and dm) we find that 



In 



H{w) _ VI + 8dw + l 
H{wa) Adw 



-In 



1+1 



/ Vl + 8dw + 1 



(1 



\-l/2 



V8d- 1 
2d 



arctan 



1 + Sdw 
8d-l 



+ .f2{wo), (42) 



where /2(ifo) is again an integration constant analogous 
to fi{wo). This solution is plotted in Fig 2] alongside the 
numerical solution. One then sees that to leading order 



H{w) 



(-1 


■ 1 


1 ■ 




\ 2d 


Wo 


w 


m 



(43) 



where the higher order terms in w. Hence, we 

have found that in the large radius Hagedorn regime the 
Hubble parameter decreases at an ever slower rate un- 
til u> ~ 1/d, beyond which it slightly increases and this 
result is independent of the initial conditions, i.e. inde- 
pendent of whether (|e|/w)initiai is smaller or greater than 
0{2d). We can further obtain the approximate time evo- 
lution of radius; it is more convenient to work in terms 
of entropy density, u; ~ a; = oq/s. After some algebraic 
manipulations ijlS]) can be integrated. The resulting in- 
complete gamma functions can be expanded about large 
s. So; keeping leading order terms only we find that 



R 

Rin 



2dao 
1 In 



2ao 



-l/d 



(44) 



The constant ao indicates that it is the momentum modes 
(as well as the oscillator modes) that play an important 
role in the evolution of the large radius universe rather 
than the winding modes. Roughly speaking, in the ab- 
sence of momentum modes, oq = 0, and hence the uni- 
verse remains static to this order. The initial entropy 
density is Sm = S/Rf^ and one can think of Rm (al- 
though not necessarily) as the radius of the universe be- 
yond which there is a negligible number of long string 
(IMI), 



Ri, 



(In 5)1/'* » 1, 



which can be thought of as the radius that distinguishes 
large from small radius Hagedorn phase. 

The characteristic timescale of the Hagedorn phase 
can be estimated from (|44p . Atn ^ e'^^" / (sinHin)- Note 
that we can further eliminate the dependence on the ini- 
tial Hubble parameter on account of (|49p in favour of the 
initial string coupling (7i„ = e*^'-. Therefore, the duration 
of the Hagedorn phase is estimated to be 



At 



1 



H 



9u 



(45) 



As one would expect, in the limit of infinite entropy 
density Sjn oo one finds that the universe is static 
R{t) = i?in while Atn oo. 

The dilaton time evolution can also be computed 
and from (gl]) and (g?]) it follows that 



In 1 



2ao 



(46) 



To see that the dilaton runs to very weak coupling notice 
that in terms of string coupling and for large times 
takes the form g ~ (7in(sin-ffinAt)~i. Then, on account 
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of p5)) and p9)) we find that at the onset of radiation 
domination the string couphng, g = e'^, wiU have reached 



gr ~ 5ine 



(47) 



which for reasonable values of the initial entropy density 
is tiny on account of (?in ^ 1- 

We shall not write down the dependence of the Ein- 
stein frame Hubble parameter on Einstein frame vol- 
ume, this amounts to substituting (|42l40l39l36p into 
It suffices to say that the Einstein frame Hubble pa- 
rameter, H, always decreases with increasing (Einstein 
frame) volume unless the dominant energy condition is 
violated, w < —1, (sec ([20]) ) or positivity of energy con- 
dition, > d, is relaxed. Furthermore, for H to in- 
crease for a large range of values for Y (rather than just 
close to y ~ —d) the condition becomes more restrictive: 
w < —{d + l)/{d— 1). These results are in accordance 
with • It is also important to note that positive energy 
and decreasing dilaton necessitate H and H to have the 
same sign, so that if the universe expands in the string 
frame it also expands in the Einstein frame. 



VI. PROBLEMS OF THE STRING GAS 

We shall now discuss various problems of the dilaton- 
gravity realisation of string gas cosmology. We shall first 
consider the issue of thermal equilibrium and improve on 
the current estimates found in [2l| (and also @) by show- 
ing that the assumption of thermal equilibrium indeed is 
not justifiable. We provide sufficient tools to actually es- 
timate the evolution of the interaction rate per Hubble 
rate as a function of size (or equivalently entropy density 
s) and present an analytic estimate of this quantity in 
the limit w 1. We then concentrate on the late time 
phenomenology of this scenario. By late time we mean 
radiation domination onwards. In what follows we make 
the assumption that the end of the Hagedorn era (where 
the universe is dominated by massless modes of string) 
is to be identified with the radiation dominated phase 
of standard Big Bang cosmology in order to extract late 
time predictions from this model. 

In the current section we show that the Hubble 
length must shrink to a tiny fraction of the size of the 
universe for a toroidal background, on which the above 
scenario is based, to be acceptable. We then go on to 
show that there is a dilaton problem associated with this 
scenario, namely that post-Hagedorn considerations re- 
quire the string coupling to be 0(1) at the onset of radia- 
tion, while the Hagedorn evolution leads to a very weakly 
coupled dilaton. In addition to these problems there is 
also an entropy problem, which stems from the fact that 
in the absence of some additional mechanism (e.g. infia- 
tion) string gas cosmology does not offer an explanation 
for the large amount of entropy in the current day uni- 
verse, although it has been speculated that the entropy 
may be generated by a period of oscillation around the 
self-dual radius (see also [H, [13, [13, H [l!] ) . 



A. Thermal Equilibrium and Jeans Instabilities 



The ques tion of thermal equilibrium has been ad- 
dressed in |2l| (see also [§|) in the deep Hagedorn regime, 
namely at the self dual point. An estimate was obtained 
and this indicated that thermal equilibrium cannot be 
maintained. However, at the self-dual point the universe 
may well be strongly coupled which leads us to distrust 
the dilaton-gravity action from which we started from. In 
particular, it is only meaningful to speak about thermal 
equilibrium at epochs where we can trust our theory. We 
shall here only consider the large volume region where 
the formalism developed above applies. The quantity of 
interest is the ratio of interaction rate per expansion rate, 
T/H. For thermal equilibrium we require this ratio to be 
much greater than 1. We can estimate the interaction 
rate at large radii from 



r 

H 



.40 



H 



(48) 



where n is the number density of strings. This was found 
in section HVl to be given by 

n{R) - l + R-'^lnS. 

We shall not display the full result but shall concentrate 
on the limiting case w <^ 1. Then, with it; ~ x = ag/s, 
with s the entropy density, the Hubble parameter (j43p 
and dilaton ((4T|) can be used to estimate the dependence 
of T/H on s{~ w/ao): 



— - Hr^e^^'" ( — 



l/2d 



-3(si„-s)/2£i 



where we have taken R > (In 3)^^'^ so that n(s) 1 as ap- 
propriate for the large radius solutions for and H. We 
can do better than this however because constrains 
the initial conditions as was emphasised in section IIIH 
^2^ ~ e2'^i"Pin/i:2^ In particular, Win < 1 and so ^ 
on account of (|39p enforces Yi^ ~ — l/win ~ — Sjn/aQ. 
Furthermore, pin = 1/xin = Sin/oo and hence the initial 
Hubble parameter must satisfy 



e'^'"(ao/sin)'/^ 



(49) 



from which we see that we can trade the initial value of 
the Hubble constant for the initial value of the dilaton. 
Substituting this into the above estimate implies that 



r 

H 



1/2 



l/2ti 



ao 



(50) 



Hence, initially, by which we mean at a radius i?in ^ 1 
(e.g. i?i„^^ (In 5)1/'^), TiJHin - e^-^-sff . So, in the 
large radius regime where the number density of long 
strings is negligible compared to the number density of 
small loops, (5 ^)i?^ > InS*, we find that the change in 
radius during which T/H ^ 1 will be 



R — Ru 

Rin 



< 



i(si„/ao) _ I 



(51) 
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We then find that to have a positive increase in radius 
with the ensemble in thermal equilibrium, T/H ^ 1, 
requires the initial string coupling to satisfy 



> 



1/6 



(52) 



the Hagedorn phase and the beginning of the radiation 
era by (and Rr in the Einstein frame), and assuming 
adiabatic expansion, we find 



Rr 



> 



To 



Ro > 10"^ mm 



(56) 



Therefore, from (|52p we learn that maintaining thermal 
equilibrium seems to require a strongly coupled initial 
configuration. 

Let us now turn to the issue of strong gravity and 
Jeans instabilities. Avoiding Jeans instabilities is prob- 
ably not a fundamental requirement: It merely must be 
satisfied to ensure that a thermodynamic treatment is 
justified. Jeans instabilities occur when gravity becomes 
strongly coupled but we do not know the law of grav- 
ity at strong coupling (e.g. there may be some non- 
perturbative strong coupling effect which prevents the 
formation of Jeans instabilities). Having said that, if we 
do require that there be no Jeans instabilities then the 
dilaton is constrained [iO] by e^'^pR^ ^ 1. We also re- 
quire i? 3> 1 for the above formalism to apply and so 
these two conditions can be satisfied provided e^'^p ^ 1. 
We can hence constrain the dilaton and the entropy den- 
sity according to 



(53) 



We are to distrust the outcome of our theory if this condi- 
tion is violated. To be specific, violation of ([55)) does not 
imply that thermodynamic equilibrium cannot be main- 
tained. This can also be viewed as an initial condition 
constraint. As we have shown above both the dilaton and 
density decrease with increasing volume so their maxi- 
mum values will in fact be their initial values. So, from 
(|5(B and (El) we find that 



H 



< s 



-3/2 



ao 



1/2 



l/2d 



-3(sin — s) /2dao 



(54) 



However, s < Si„ and so F/iJ <C ^ 1. Therefore, in 
the domain of validity of the above formalism the string 
ensemble cannot be in thermal equilibrium; we have how- 
ever seen some hints that an initially strongly coupled 
dilaton may have an important role to play. 



where and Tq are the temperatures at the start of 
the radiation era and today respectively, with Tr < mpi- 
Hence the universe at the end of the Hagedorn era is very 
large in string units, justifying our concentrating on the 
large-radius thermodynamics. 

The fact that the universe is large can be re- 
expressed as problem for the dilaton, as follows. Working 
in the Einstein frame and taking the number of large di- 
mensions d = 3, we can estimate the ratio, RrHr, at 
time tr, when the FRW radiation era (constant dilaton) 
begins, given that 




where i?cq is the radius at matter-radiation equality. 
Note also that we are neglecting the dark energy era for 
our rough estimate. It then follows that 



RrHr — RqHq 



Rq 



Given Rcq/Ro ~ 10 and that the entropy density s oc 
R^^ in the FRW cosmology, we have 

RrHr > 10-2 10-2 ( 




So 



Tn 



Therefore, at the end of the Hagedorn regime this ratio 
should satisfy 



RrHr — RrHr ^ 10 



29 / '^r 



(57) 



B. The Size Problem 

At present we know from measurements of the CMB 
power spectrum at large angular scales that if the Uni- 
verse is toroidal, its radius, Rq, cannot be much smaller 
than the Hubble length today, H^^, (see e.g. [4l|): 



RqHq = RqHq > 1, 



(55) 



where we have used the fact that the string frame ratio 
RH is equal to the Einstein frame ratio RH when the 
dilaton is constant. Denoting the radius at the end of 



assuming Tr ^ ?Hag — i^pi- Hence, at the beginning of 
the radiation era, independently of whether we are work- 
ing in the string or Einstein frame, the Hubble length 
must occupy a tiny fraction of the size of the Universe 
for a toroidal background to be phenomenologically ac- 
ceptable. 

We can estimate the ratio on account of the 

string frame expression (|43p with w ~ ao/s, Sh = sH^'^ 
and S = sR"^, given that at the onset of radiation s 1, 



RrHr ^ 4^''^ exp 



2dac 



RinHi, 



(58) 
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Therefore, one way of satisfying the constraint ((57)) is for 
the initial (string frame) Hubble parameter, 




l,-l/2gSi„/(6ao) 



i?in > 10 



where we have taken d — S. For temperatures ~ mpi 
and given that Sm 3> 1, this is much larger than that 
allowed from the assumption of adiabaticity according to 
which < 1. On account of the initial conditions con- 
straint ([49]) one sees that this inequality can be rewritten 
in terms of the initial string coupling, g^, 



9in > 10 



29 



JTlpi 



p--lg.s.„/(6ao) 



which is inconsistent with requirement of weak coupling. 
We conclude that the 3-torus (on which our universe is 
embedded in this scenario) cannot become large enough 
to agree with observations, namely equation ()55p . if 
string is weakly coupled close to the string scale, gin <C 1. 



C. The Dilaton Problem 

A dilaton problem arises in the following way. We 
have seen that the onset of radiation domination is char- 
acterised by a string frame entropy density of order one, 
Sr = S/R^ ^ 1. We can derive the entropy density in 
the Einstein frame, s = S/R'^, on account oi R ~ g^^R 
which relates Einstein and string frame radii ([T8|). We 
then find that the Einstein frame entropy density at the 
onset of radiation is 

Sr ~ gl- (59) 
The entropy density is approximately s ^ and hence 

9r ^ , 

mpi 

where we have explicitly restored the units. To a first 
approximation one expects the temperature at the end 
of the Hagcdorn region to be Tr ^ O(mpi) and hence we 
have shown that the string coupling should be of order 
one, 



gr ^ 1. 



(60) 



This we expect independently from the relation between 
gauge and string couplings in a universe where the com- 
pact dimensions remain at the string scale. However, as 
we have seen above, during the Hagedorn era the dilaton 
is expected to run to negligible coupling and at the onset 
of radiation (l47l). 



gr - ginS ""^ < 1, (61) 

where g^n ^1- It is hence clear that there is large dis- 
agreement between the string coupling at the onset of ra- 
diation domination as required from the post-Hagedorn 
considerations (j60p and that predicted from the equa- 
tions of motion during the large radius Hagedorn regime 

m. 



VII. DISCUSSION 

We have considered an ensemble of high energy 
strings in dilaton gravity background. The picture we 
have is the following: the universe is initially at string 
scale, all dimensions are of equal size and all degrees of 
freedom are present. Quantum fluctuations create a hier- 
archy of scales and based on the original formulation [l[ 3 
dimensions will become large. This is our starting point 
for the calculations presented in this paper, 3 dimensions 
are large and 6 remain at string scale. 

We have calculated the equation of state for a (het- 
erotic) string gas and have used it to determine the result- 
ing dynamics of the metric and dilaton fields. We have 
concentrated on the large radius Hagedorn regime with 
totally compact topology and have allowed d of the 9 — d 
spatial dimensions to freely evolve. We have found ana- 
lytic solutions for both the Hubble parameter and dilaton 
as functions of size for the entire large radius Hagedorn 
evolution and have compared our results with numerical 
solutions. We have found that the equation of state has 
the form 

w{x) 



1-x' 

with x/aQ the inverse entropy density, s, and oq a con- 
stant arising due to the presence of momentum modes. 
Higher order terms can be found in (j25p . For large en- 
tropy densities, s, we have found that (j){s) ~ 0in — 
(sin - s)/2dao, and H{s) ~ i?in(s/sin)^/^'' exp{-(si„ - 
s)/2dao}, the full expressions arc given in ([3D]) and ([32]) 
respectively. In particular, the dilaton is shown to de- 
crease inversely proportional to the volume of the large 
compact dimensions. The Hubble parameter decreases 
with increasing volume, reaches a minimum when the 
equation of state w ~ 1/d and then may increase slightly, 
this is sensitive to the dependence of the total density of 
states, n, on E and R when the energy density p is close 
to the Hagedorn energy density pi . Integrating the lead- 
ing order form for the Hubble parameter as a function of 
size we have derived the leading dependence of size with 
time ([33]) and have obtained an estimate for the duration 
of the Hagedorn era in the large radius Hagedorn regime 



At 



1 



H 



gin Y -^in 



which is very large if the dilaton is weakly coupled. 

It is also useful to compare our results for radius 
and dilaton as functions of time for large radius evolu- 
tion with the corresponding results for vanishing pres- 
sure at small radius. A main difference between 
small and large radius evolution is that at small radii 
the pressure is very small and slowly evolving while at 
large radii the pressure is approximately constant and of 
0(1). It suffices to say that at large times TV Q find 
that (j){t) ~ ~\nt + const. Our solution ([ID]) agrees in 
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that we have also found a logarithmic dependence of dila- 
ton on time, (j>{t) ~ — ln(l + giy^y/'si^At/2) + const. As 
far as the radius is concerned we have found that there 
is no maximum radius due to the fact that the equa- 
tions of state increases with increasing size. In partic- 
ular, TV [2] find that A = Ai„ + ln(f^)i/^ for some 
integration constant c whereas we find that in the large 
radius region A ~ Ai„ + ln[l-7ln(l + i7i„Ai/7)]~^/'' with 
7 = 2dao/sni <C 1 vahd for At < e*'"/(5in\/sin) (use has 
also been made of (|49l) ). 

Therefore, the results we have obtained differ from 
the small radius evolution. The source of these differ- 
ences boil down to the string spectrum. In particular, 
even though the oscillator modes account for most of the 
energy of the ensemble, the number of long strings at 
large radii, S^^'^ 3> i? S> {InSY^'^, is neghgible com- 
pared to the number of small loops present. The two 
types of string contribute equally to the number density 
when the radius R ~ (InS')^/''. Furthermore, the mass- 
less modes are expected to dominate the ensemble below 
the Hagedorn energy density [l| where R ~ S^/'^. If 
above the radius R ~ (In Sy^'^ winding strings are few in 
number then one would naively expect that the splitting 
from 9 to 6-1-3 dimensions (in the original manifestation 
of the Brandenberger-Vafa mechanism of spacetime di- 
mensionality) should occur before this radius is reached. 



However, as we have stressed above evolution in the small 
radius regime is very speculative to date. The main rea- 
son being that one expects there to be non-perturbative 
effects that come into play and these may in turn give 
rise to a vastly different picture from the one obtained 
by perturbative calculations alone. 

In the last section we have seen that there are a 
number of problems associated with a large radius high- 
entropy string gas universe, associated with the lack of 
thermodynamic equiilibrium and the requirement to have 
a large but not too weakly-coupled universe today. One 
cannot rely on early strongly-coupled dynamics to re- 
trieve the scenario, as a small radius universe must be 
followed by a large radius one, with large entropy if we 
are to describe our universe, and our analysis applies to 
this phase. This version of string gas cosmology therefore 
still has many difficulties to overcome. 
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